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Abstract 

The unicellular green algae Chlamydomonas swims with two flagella, which can synchronize their beat. 
Synchronized beating is required to swim both fast and straight. A long-standing hypothesis proposes 
that synchronization of flagella results from hydrodynamic coupling, but the details are not understood. 
Here, we present realistic hydrodynamic computations and high-speed tracking experiments of swimming 
cells that show how a perturbation from the synchronized state causes rotational motion of the cell body. 
This rotation feeds back on the flagellar dynamics via hydrodynamic friction forces and rapidly restores 
the synchronized state in our theory. We calculate that this 'cell body rocking' provides the dominant 
contribution to synchronization in swimming cells, whereas direct hydrodynamic interactions between the 
flagella contribute negligibly. We experimentally confirmed the coupling between flagellar beating and cell 
body rocking predicted by our theory. We propose that the interplay of flagellar beating and hydrodynamic 
forces governs swimming and synchronization in Chlamydomonas. 
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Eukaryotic cilia and flagella are long, slender cell appendages that can bend rhythmically and 
thus represent a prime example of a biological oscillator [ 1 ] . The flagellar beat is driven by the col- 
lective action of dynein molecular motors, which are distributed along the length of the flagellum. 
The beat of flagella, with typical frequencies ranging from 20 — 60 Hz, pumps fluids, e.g. mucus 
in mammalian airways [2], and propels unicellular micro- swimmers like Paramecium, spermato- 
zoa, or algae [3]. The coordinated beating of collections of flagella is important for efficient fluid 
transport [2, 4] and fast swimming [5]. This coordinated beating represents a striking example 
for the synchronization of oscillators, prompting the question of how flagella couple their beat. 
Identifying the specific mechanism of synchronization can be difficult as synchronization may oc- 
cur even for weak coupling [6]. Further, the effect of the coupling is difficult to detect once the 
synchronized state has been reached. 

Hydrodynamic forces were suggested to play a significant role for flagellar synchronization 
already in 1951 by G.I. Taylor [7]. Since then, direct hydrodynamic interactions between flagella 
were studied theoretically as a possible mechanism for flagellar synchronization [8-11]. Another 
synchronization mechanism that is independent of hydrodynamic interactions was recently de- 
scribed in the context of a minimal model swimmer [12]. This mechanism crucially relies on the 
interplay of swimming motion and flagellar beating. 

Here, we address the hydrodynamic coupling between the two flagella in a model organism 
for flagellar coordination [13-16], the unicellular green algae Chlamydomonas. Chlamydomonas 
propels its ellipsoidal cell body, which has typical diameter 10 fxm, using a pair of flagella, whose 
lengths are about 10 jum [13]. The two flagella beat approximately in a common plane, which is 
collinear with the long axis of the cell body. In that plane, the two beat patterns are nearly mirror- 
symmetric with respect to this long axis. The beating of the two flagella of Chlamydomonas can 
synchronize, i.e. adopt a common beat frequency and a fixed phase relationship [13-16]. bi- 
phase synchronization of the two flagella is required for swimming along a straight path [16]. The 
specific mechanism leading to flagellar synchrony is unclear. 

Here, we use a combination of realistic hydrodynamic computations and high-speed track- 
ing experiments to reveal the nature of the hydrodynamic coupling between the two flagella of 
free swimming Chlamydomonas cells. Previous hydrodynamic computations for Chlamydomonas 
used either resistive force theory [17, 18], which does not account for hydrodynamic interactions 
between the two flagella, or computationally intensive finite element methods [19]. We employ 
an alternative approach and represent the geometry of a Chlamydomonas cell by spherical shape 
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primitives, which provides a computationally convenient method that fully accounts for hydrody- 
namic interactions between different parts of the cell. Our theory characterizes flagellar swimming 
and synchronization by a minimal set of effective degrees of freedom. The corresponding equa- 
tion of motion follows naturally from the framework of Lagrangian mechanics, which was used 
previously to describe synchronization in a minimal model swimmer [12]. These equations of 
motion embody the basic assumption that the flagellar beat speeds up or slows down according to 
the hydrodynamic friction forces acting on the flagellum. This assumption is supported by pre- 
vious experiments, which showed that the flagellar beat frequency depends on the viscosity of 
the surrounding fluid if the concentration of the chemical fuel ATP is sufficiently high [20]. The 
simple force- velocity relationship for the flagellar beat employed by us coarse-grains the behavior 
of thousands of dynein molecular motors that collectively drive the beat. Similar force-velocity 
properties have been described for individual molecular motors [21] and reflect a typical behavior 
of active force generating systems. 

Our theory predicts that any perturbation of synchronized beating results in a significant yaw- 
ing motion of the cell, reminiscent of rocking of the cell body. This rotational motion imparts 
different hydrodynamic forces on the two flagella, causing one of them to beat faster and the other 
to slow down. This interplay between flagellar beating and cell body rocking rapidly restores 
flagellar synchrony after a perturbation. Using the framework provided by our theory, we analyze 
high-speed tracking experiments of swimming cells, confirming the proposed two-way coupling 
between flagellar beating and cell body rocking. 

Previous experiments restrained Chlamydomonas cells from swimming, holding their cell body 
in a micropipette [14-16]. Remarkably, flagellar synchronization was observed also for these con- 
strained cells. This observation seems to argue against a synchronization mechanism that relies 
on swimming motion. However, the rate of synchronization observed in these experiments was 
faster by an order of magnitude than the rate we predict for synchronization by direct hydrody- 
namic interactions between the two flagella in the absence of any motion. In contrast, we show 
that rotational motion with a small amplitude of a few degrees only, which may result from either 
a residual rotational compliance of the clamped cell or an elastic anchorage of the flagellar pair, 
provides a possible mechanism for rapid synchronization. 
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I. RESULTS AND DISCUSSION 



A. High-precision tracking of confined Chlamydomonas cells 

To study the interplay of flagellar beating and swimming motion, we recorded single wild-type 
Chlamydomonas reinhardtii cells swimming in a shallow observation chamber using high-speed 
phase-contrast microscopy (1000 fps). The chamber heights were only slightly larger than the cell 
diameter so that the cells did not roll around their long body axis, but only translated and rotated 
in the focal plane. This confinement of cell motion to two space dimensions and the fact that 
the approximately planar flagellar beat was parallel to the plane of observation greatly facilitated 
data acquisition and analysis. From high-speed recordings, we obtained the projected position and 
orientation of the cell body as well as the shape of the two flagella, see figure 1A. 

In a reference frame of the cell body, each flagellum undergoes periodic shape changes. To for- 
malize this observation, we defined a flagellar phase variable by binning flagellar shapes according 
to shape similarity, see figure IB. A time-series of flagellar shapes is represented by a point cloud 
in an abstract shape space. This point cloud comprises an effectively one-dimensional shape cy- 
cle, which reflects the periodicity of the flagellar beat. Each shape point can be projected on the 
centerline of the point cloud. We define a phase variable (p running from to 2tt that parametrizes 
this limit cycle by requiring that the phase speed ifi is constant for synchronized beating. Approx- 
imately, we determine this parametrization from the condition that the averaged phase speed is 
independent of the location along the limit cycle. This defines a unique flagellar phase for each 
tracked flagellar shape. The width of the point cloud shown in figure IB is a measure for the vari- 
ability of the flagellar beat during subsequent beat cycles. We find that the variations of flagellar 
shapes for the same value of the phase variable are much smaller than the shape changes during 
one beat cycle. For our analysis, we therefore neglect these variations of the flagellar beat. In 
this way, we characterize a swimming Chlamydomonas cell by 5 degrees of freedom: its position 
(x, y) in the plane, the orientation angle a of its cell body, and the two flagellar phase variables ip L 
and ifR for the left and right flagellum, respectively. Our theoretical description will employ the 
same 5 degrees of freedom and use flagellar shapes tracked from experiment for the hydrodynamic 
computations. 
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B. Hydrodynamic forces and interactions 

For a swimming Chlamydomonas cell, inertial forces are negligible (as characterized by a low 
Reynolds number of Re ~ 1CT 3 [19]), which implies that the hydrodynamic friction forces exerted 
by the cell depend only on its instantaneous motion [22] . To conveniently compute hydrodynamic 
friction forces and hydrodynamic interactions, we represented the geometry of a Chlamydomonas 
cell by N = 300 spherical shape primitives, see figure 2A. The spheres constituting the cell body 
are treated as a rigid cluster. For simplicity, we consider free swimming cells and do not include 
wall effects in our hydrodynamic computations. Flagellar beating and swimming corresponds to 
a simultaneous motion of all 300 spheres of our cell model. The dependence of the corresponding 
hydrodynamic friction forces and torques on the velocities of the individual spheres is charac- 
terized by a grand hydrodynamic friction matrix G. We computed this friction matrix G using 
a Cartesian multipole expansion technique [23], see 'Materials and Methods' for details. Fig- 
ure 2C,D shows a sub-matrix that relates force and velocity components parallel to the long axis 
of the cell. The entries of the color matrix depict the force exerted by any of the flagellar spheres 
or by the cell body cluster (row index), if a single flagellar sphere or the cell body cluster is moved 
(column index). The indexing of flagellar spheres is indicated by cartoon drawings of the cell next 
to the color matrix. The diagonal entries of this friction matrix are positive and account for the 
usual Stokes friction of a single 'flagellar sphere' (or of the cell body). Off-diagonal entries are 
negative and represent hydrodynamic interactions. We find considerable hydrodynamic interac- 
tions between spheres of the same flagellum, as well as between each flagellum and the cell body. 
However, interactions between the two flagella are comparably weak. 

C. Theoretical description of flagellar beating and swimming 

We present dynamical equations for the minimal set of 5 degrees of freedom shown in figure 1 A 
to describe flagellar beating, swimming, and later flagellar synchronization in Chlamydomonas. 
These equations of motion follow naturally from the framework of Lagrangian mechanics of dis- 
sipative systems, which defines generalized forces conjugate to effective degrees of freedom. 

Motivated by our experiments, we describe the progression through subsequent beat cycles of 
each of the two flagella by respective phase angles (p L and ip R , see figure 1A. The angular fre- 
quency Uj of flagellar beating is given by the time-averaged phase speed so we can think of 
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the phase speed as the instantaneous beat frequency. We are interested in variations of the phase 
speed that can restore a synchronized state after a perturbation. We introduce the key assumption 
that changes in hydrodynamic friction during the flagellar beat cycle can increase or decrease the 
phase speed of each flagellum. Specifically, we assume that for both the left and right flagellum, 
j = L, R, the respective flagellar phase speed ifij is determined by a balance of an active driving 
force Qj that coarse-grains the active processes within the flagellum and a generalized hydrody- 
namic friction force Pj, which depends on ipj. Note that in addition to hydrodynamic friction, 
dissipative processes within the flagella may contribute to the friction forces P L and P R . We do 
not consider such internal friction in our description as it does not change our results qualitatively. 
The hydrodynamic friction forces Pj have to be computed self-consistently for a swimming cell. 
We restrict our analysis to planar motion in the 'xy' -plane and thus consider the position (x, y) 
and the orientation a of the cell body with respect to a fixed laboratory frame, see figure 1A. 

Any change of the degrees of freedom x, y, a, (p L , (p R results in the dissipation of energy into 
the fluid at some rate 1Z. This dissipation rate 1Z characterizes the mechanical power output of 
the cell and plays the role of a Rayleigh dissipation function known in Lagrangian mechanics; 
it can be written as 1Z = xP x + yP y + dP a + ^lPl + ^rPr, which defines the generalized 
friction forces conjugate to the different degrees of freedom. The forces P L , P R , P a are conjugate 
to an angle and have physical unit pN /im. We compute the generalized friction forces using the 
grand hydrodynamic friction matrix G introduced above. In brief, the superposition principle of 
low Reynolds number hydrodynamics relevant for Chlamydomonas swimming [22] implies that 
the generalized friction forces relate linearly to the generalized velocities, Pj = T jx x + T jy y + 
Tj a a+Tj L (p L +Tj R (p R . This defines the generalized hydrodynamic friction coefficients T^, which 
are suitable linear combinations of the entries of the grand hydrodynamic friction matrix G, see 
'Materials and Methods'. 

The friction force P x conjugate to the re-coordinate of the cell position represents just the x- 
component of the total force exerted by the cell on the fluid, an analogous statement applies for P y ; 
P a is the total torque associated with rotations around an axis normal to the plane of swimming. If 
the swimmer is free from external forces and torques, we have P x = P y = and P a = 0. Together 
with the proposed balance of flagellar friction and driving forces, Pl = Ql and P R = Q R , we 
have a total of 5 force balance equations, which allow us to solve for the time derivatives of the 5 
degrees of freedom. We obtain an equation of motion that combines swimming and flagellar phase 
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dynamics 

(x, y, a, <p L , <p R ) T = T-\0, 0, 0, Q L , Q R ) T . (1) 

The phase dependence of the active driving forces Qj(<fij) is uniquely specified by the condition 
that the phase speeds should be constant, (fij = u , for synchronized flagellar beating with zero 
flagellar phase difference 5 = 0, where 5 = <£l — Vr- 

In essence, this generic description implies that the phase speed of one flagellum is determined 
by hydrodynamic friction forces, which in turn depend on the swimming motion of the cell. Since 
the swimming motion is determined by the beating of both flagella, equation 1 effectively defines 
a feedback loop that couples the two flagellar oscillators. 

D. Theory and experiment of Chlamydomonas swimming 

Using the equation of motion (equation 1), we can compute the swimming motion of our model 
cell. For mirror- symmetric flagellar beating with zero flagellar phase difference 5 = 0, the model 
cell follows a straight path with an instantaneous velocity that is positive during the effective 
stroke, but becomes negative during a short period of the recovery stroke (figure 3 A, left panel). 
The cell swims two steps forward, one step back. This saltatory motion is also observed exper- 
imentally (figure 3 A, middle panel). In our computation, the instantaneous swimming velocity 
reaches values up to 200 /im/s, which agrees with experimental measurements for free-swimming 
cells [24], but overestimates the observed translational swimming speeds in shallow chambers, in 
which wall effects are expected to reduce the speed of translational motion (compare left and right 
panels in figure 3A). If the two flagella are beating out of phase, the cell will not swim straight 
anymore, but the cell body yaws, see figure 3B. Cell body yawing is observed experimentally 
(right panel), with measured yawing rates that agree well with our computations (left panel). The 
proximity of boundary walls is known to reduce translational motion, but to affect rotational mo- 
tion to a much lesser extent for a given distance from the wall [18]. This is indeed observed in our 
experiments with cells swimming in shallow chambers: while the observed translational speed is 
smaller than predicted (figure 3A), the observed yawing rates are very similar to the predicted ones 
(figure 3B). The good agreement between theory and experiment for the yawing rate supports our 
hydrodynamic computation as well as our description of flagellar beating using a single phase vari- 
able. In the next section, we show that rotational motion is crucial for flagellar synchronization, 
while translational motion is less important. 
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E. Theory of flagellar synchronization by cell body yawing 

We now demonstrate how yawing of the cell body leads to flagellar synchronization. We first 
examine the flagellar phase dynamics after a perturbation of in-phase flagellar synchrony. Figure 
4A shows numerical results for a free swimming cell obtained from solving the equation of motion 
(equation 1). The initial flagellar asynchrony causes a yawing motion of the model cell, which is 
characterized by periodic changes of the cell's orientation angle a(t). The phase difference 5 
between the left and right flagellum decays approximately exponentially as 5(t) ~ exp(— Xt/T) 
with a rate constant A [measured in beat periods T = 2ir/u ] that will serve as a measure of the 
strength of synchronization. 

To mimic experiments where external forces constrain cell motion, we now consider the ide- 
alized case of a cell that cannot translate, while cell body yawing is constrained by an elastic 
restoring force Q a = — ka. Again, the two flagella synchronize in-phase, provided some residual 
cell body yawing is allowed, see figure 4B. In the absence of an elastic restoring force (k = 0), 
when the model cell cannot translate, but can still freely rotate, its yawing motion and synchro- 
nization behavior is very similar to the case of a free swimming cell that can rotate and translate. 
For a fully clamped cell body, however, the synchronization strength is strongly attenuated, and 
is solely due to the direct hydrodynamic interactions between the two flagella. In this case of 
synchronization by hydrodynamic interactions, the time-constant for synchronization is decreased 
approximately 20-fold compared to the case of free swimming. These numerical observations 
point at a crucial role of cell body yawing for flagellar synchronization. The underlying mecha- 
nism of synchronization can be explained as follows. For in-phase synchronization, the flagellar 
beat is mirror- symmetric and the cell swims along a straight path. If, however, the left flagellum 
has a small head-start during the effective stroke, this causes a counter-clockwise rotation of the 
cell, see figure 3B. This cell body yawing increases (decreases) the hydrodynamic friction encoun- 
tered by the left (right) flagellum, causing the left flagellum to beat slower and the right one to beat 
faster. As a result, flagellar synchrony is restored. 

Next, we present a formalized version of this argument using a reduced equation of motion. We 
thus arrive at a simple theory for biflagellar synchronization, which will later allow for quantitative 
comparison with experiments. As in figure 4B, we assume that the cell is constrained such that 
it cannot translate (x — y — 0). The cell can still yaw, possibly being subject to an elastic 
restoring force Q a = —ka. This leaves only three degrees of freedom: (p L , (p R , and a. Neglecting 
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direct hydrodynamic interactions between the flagella, we can reduce the equations of motion for 
a clamped cell (equation 1 with constraint x = y = 0) to a set of three coupled equations for the 
three remaining degrees of freedom 

<-Pl = oj l - ii((fL)a, (2) 
p R = uj R + p(p R )a, (3) 
k a + p(ip L , ip R ) a = -v(p L )p L + v(<p R )<p R . (4) 

The coupling function p in equation 2 characterizes the effect of cell body yawing on the flagellar 
beat as detailed below, while v describes how asynchronous flagellar beating results in yawing; p 
is the hydrodynamic friction coefficient for yawing of the whole cell. The coupling functions p, 
v, and p can be computed using our hydrodynamic model[35]. Their dependence on the flagellar 
phase is shown in figure 5 (left panels). The physical significance of equations 2-4 can be explained 
as follows: Equation 2 implies that during the effective stroke of the left flagellum (ip ~ 0°), a 
counter-clockwise rotation of the whole cell slows down the flagellar beat, while a clockwise 
rotation speeds it up (figure 5B, p > 0). Equation 3 implies the converse for the right flagellum. 
During the recovery stroke (p ~ 180°), the effect is opposite and a counter-clockwise rotation 
of the cell would speed up the beat of the left flagellum (p < 0). Equation 4 states that flagellar 
beating causes the cell body to yaw: if the right flagellum were absent, the model cell rotates 
clockwise (a < 0) during the effective stroke of the left flagellum (figure 5 A, v > 0), and counter- 
clockwise during its recovery stroke {y < 0). This swimming behavior is observed for uniflageilar 
mutants [18]. For synchronized beating of the two flagella, the right-hand side of equation 4 
cancels to zero and the model cell swims straight. For asynchronous flagellar beating with a finite 
phase-difference 5 = p>L — ^r, the phase-dependence of the coupling function u(p) results in an 
imbalance of the torques generated by the left and right flagellum, respectively, which is balanced 
by a rotation of the whole cell. A similar situation arises for an elastically anchored flagellar 
pair as detailed in the Supporting Information (SI). In that scenario, the flagellar pair pivots as a 
whole when the two flagella beat asynchronously, which causes rapid synchronization analogous 
to synchronization by cell body yawing considered here. 

We study the dynamical system given by equations 2-4 after a small perturbation of the syn- 
chronized state at t = with initial flagellar phase difference < 5(0) «C 1. For simplicity, we 
assume intrinsic beat frequencies, uo L = u R = u . For small initial perturbations, the flagellar 
phase difference will either decay or grow exponentially, if we average its dynamics over a beat 
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cycle, S(t) 5(0) exp(— Xt/T). A positive value of the synchronization strength A implies that 
the perturbation decays, and that the in-phase synchronized state with 5 = is stable. The syn- 



where a prime denotes differentiation with respect to (p. Using the coupling functions p, v, p 
computed above, we obtain A > 0, which implies stable in-phase synchronization, see figure 4. 

In the case of a stiff elastic constraint, we obtain a different result for the synchronization 
strength A 



Synchronization in the absence of an elastic restoring force as characterized by equation 5, 
and synchronization involving a strong elastic coupling as characterized by equation 6 show in- 
teresting differences, which relate to the fact that in the first case the flagellar phase dynamics 
depends only on the yawing rate a, but not on a itself. The difference between these two syn- 
chronization mechanisms is best illustrated in a special case, in which both the ratio a — vjp and 
p are constant. A constant a correspond to an active flagellar driving force that does not depend 
on the flagellar phase, whereas for constant p, the angular friction for yawing would not depend 
on the flagellar configuration. In the limit of a stiff elastic constraint, k 3> pcuo, we readily find 
A = —acu <fisis"/k = auj §(v') 2 jk > 0, which indicates stable in-phase synchronization. In 
the limit of a weak elastic constraint, k <C pu , however, the integral on the right-hand side of 
equation 5 evaluates to zero, which implies that synchronization does not occur. Hence, synchro- 
nization in the absence of an elastic restoring force requires that either vjp or p depend on the 
flagellar phase. 

For our realistic Chlamydomonas model, p and v differ (figure 5A), and also p is not constant 
(figure S3 in the SI). This allows for rapid synchronization also in the absence of elastic forces. 
Previous work on synchronization in minimal systems showed that elastic restoring forces can fa- 
cilitate synchronization [10, 25]. Here, we have shown that elastic forces can increase the synchro- 
nization strength (figure 4), but they are not required for flagellar synchronization in swimming 
Chlamydomonas cells, even if hydrodynamic interactions are neglected. 

Our discussion of flagellar synchronization can be extended to the case, where the intrinsic 
beat frequencies of the two flagella do not match. If the frequency mismatch \cu L — cu R \ is small 



chronization strength is given by A = — J Q T dt8/8. In the limit of a small elastic constraint, we 
find (see SI for details) 




(5) 




(6) 
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compared to the inverse time-scale of synchronization A/T, a general result implies that the two 
flagellar oscillators will still synchronize [6]. For a frequency mismatch that is too large, however, 
the two flagella will display phase drift with monotonously increasing phase difference. 

F. Experiments show coupling of beating and yawing 

We reconstructed the coupling functions fi(ip) and v(<p) between beating and yawing from ex- 
perimental data using the theoretical framework developed in the previous section. In brief, (1) we 
extracted the instantaneous yawing rate a and flagellar phase speeds L and ip R from high-speed 
videos of swimming Chlamydomonas cells, (2) we represented the coupling functions by a trun- 
cated Fourier series, and (3) we obtained the unknown Fourier coefficients by linear regression 
using equations 2-4. The high temporal resolution of our imaging enabled us to accurately deter- 
mine phase speeds as time-derivatives of flagellar phase angle data. Figure 5B displays averaged 
coupling functions obtained by fitting for a typical Chlamydomonas cell. We find a significant cou- 
pling between flagellar phase speeds and yawing rates, which are in good qualitative agreement 
with the theoretical predictions. Figure S7 in the SI shows fits for 5 more cells. 

For the experimental conditions used, we commonly observed cells that displayed a large fre- 
quency mismatch between the two flagella. For the cells selected for analysis, this frequency 
mismatch exceeded 30%. This large frequency mismatch caused flagellar phase drift, making 
these cells ideally suited to study flagellar coupling. Experiments were done using either white 
light illumination, which gave maximal image quality, or red light illumination, which reduces 
a possible phototactic stimulation of the cells. In all cells, we consistently observed a coupling 
of flagellar beating and cell body yawing for both red and white illumination (figure S7). This 
coupling should be a manifestation of the internal working mechanics of the flagellum and the 
hydrodynamics of swimming, and hence represent a general feature of mechanical regulation of 
flagellar dynamics. We therefore argue that this coupling, which we observed in cells that do not 
synchronize, is responsible for synchronization in cells with smaller mismatch of beat frequencies. 

H. CONCLUSION AND OUTLOOK 

We have presented a theory on the hydrodynamic coupling underlying flagellar synchroniza- 
tion in swimming Chlamydomonas cells. We have shown that direct hydrodynamic interactions 



11 



between the two flagella as considered in refs. [8-10] give only a minor contribution to the com- 
puted synchronization strength and are unlikely to account for the rapid synchronization observed 
in experiments [13-16]. In contrast, rotational motion of the swimmer caused by asynchronous 
beating imparts different hydrodynamic friction forces on the two flagella, which rapidly brings 
them back in tune: Chlamydomonas rocks to get into synchrony. 

Using high-speed tracking experiments, we could confirm the two-way coupling between flag- 
ellar beating and cell body yawing predicted by our theory. The striking reproducibility of our 
fits for the corresponding coupling functions and their favorable comparison to our theory is 
highly suggestive for a regulation of flagellar phase speed by hydrodynamic friction forces that 
depends on rotational motion. Thus, coupling of flagellar beating and cell body yawing provides a 
strong candidate for the mechanism that underlies flagellar synchronization of swimming Chlamy- 
domonas cells. 

To explain a previously observed synchronization for cells held in a micropipette [14-16], we 
propose a finite clamping compliance that still allows for residual cell body yawing with an am- 
plitude of a few degrees, which is sufficient for rapid synchronization. Alternatively, a compliant 
basal anchorage of the flagellar pair or bending deformations of the elastic cell body would allow 
for flagellar synchronization by a completely analogous mechanism. The simple theory for biflag- 
ellar synchronization by rotational motion presented in this manuscript (equations 2-4) is generic 
and applies also to the pivoting motion of an elastically anchored flagellar pair as shown in the 
Supporting Information. From the observed value A = 0.3 for the synchronization strength in 
clamped cells [16], we estimate a rotational stiffness of k ~ 10 4 pN/xm for either of these two 
cases. 

Finally, the coupling of two phase oscillators by a third degree of freedom, in this case rotational 
motion, could allow for synchronization also in other contexts. For example, one may consider 
that synchronization in ciliar arrays [2] is mediated by an elastic coupling through the matrix with 
elastic deformations playing the role of the third degree of freedom. 

III. HYDRODYNAMIC COMPUTATION OF SWIMMING CHLAMYDOMONAS 

We represent a Chlamydomonas cell by an ensemble of 300 spheres of radius a = 0.25 fira, see 
Figure 2A, and use a freely available hydrodynamic library based on a Cartesian multipole expan- 
sion technique [23] to compute the grand hydrodynamic friction matrix G [22] for this ensemble 



12 



of spheres. We assume a rigid cell body, and hence that the spheres constituting the cell body 
move as a rigid unit, which results in n = 2 • 14 + 1 independently moving objects. The matrix G 
has dimensions 6n x 6n and relates the components of the translational and rotational velocities, 
Vj and Cli, of each of the n objects to the hydrodynamic friction forces and torques, Fj and Tj, ex- 
erted by the j-th object on the fluid, (F lx , F ly , F lz ,T lx , T ly , T lz , F 2x , F 2y , . . . , T ny , T nz ) = Gq 
with q = {v\ x , vi y , viz, £lix, &iy, &iz, v 2x , v 2y , ■ ■ ■ , £l ny , Qnz)- Figure 2C,D shows the sub-matrix 
Gij :V y = G ei - 4 £j-4, which relates ^-components of velocities and y-components of hydrody- 
namic forces. The reduced friction matrix T for a set of m effective degrees of freedom q is 
computed from G as T = L T GL with Qn x m transformation matrix = dq ,i/dqj, where 
q = (x, y, a,(fL, ¥>r) [12]. Initial tests confirmed that the friction matrix of only the cell body 
gave practically the same result as the analytic solution for the enveloping spheroid; similarly the 
computed friction matrix of only a single flagellum matched the prediction of resistive force theory 
[22]. 

IV. IMAGING CHLAMYDOMONAS SWIMMING IN A SHALLOW OBSERVATION CHAMBER 

Cell Culture: Chlamydomonas reinhardtii cells (CC-125 wild type mt+ 137c, R.P. Levine via 
N.W. Gillham, 1968) were grown in 300 ml TAP+P buffer [26] (with 4 x phosphate) at 24° for 2 
days under conditions of constant illumination (2x75 W, fluorescent bulb) and constant air bub- 
bling to a final density of 10 6 cells/ml. 

High Speed Video Microscopy: An assay chamber was made of pre-cleaned glass and sealed 
using Valap, a 1:1:1 mixture of lanolin, paraffin, and petroleum jelly, heated to 70° C. The surface 
of that chamber was blocked using casein solution (solution of casein from bovine milk 2 mg/ml, 
for 10 min) prior to the experiment. Single, non interacting cells were visualized using phase 
contrast microscopy set up on a Zeiss Axiovert 100 TV Microscope using a 63x Plan-Apochromat 
NA1.4 PH3 oil lens in combination with an 1.6x tube lens and an oil phase-contrast condenser 
NA 1.4. The sample was illuminated using a 100 W Tungsten lamp. For red-light imaging, an e- 
beam driven luminescent light pipe (Lumencor, Beaverton, USA) with spectral range of 640-657 
nm and power of 75 mW was used. The sample temperature was kept constant at 24° C using an 
objective heater (Chromaphor, Oberhausen, Germany). For image acquisition, an EoSens Cmos 
high-speed camera was used. Videos were acquired at a rate of 1000 fps with exposure times of 
1 ms (white light) and 0.6 ms (red light). Finally, cell positions and flagellar shapes were tracked 
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using custom-build Matlab software, see SI for details. 
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V. FIGURES 




FIG. 1: Five degrees of freedom for Chlamydomonas. A. In our experiments, conducted in shallow obser- 
vation chambers, Chlamydomonas cells swim in a plane. At each time, the position and orientation of the 
cell body is characterized by its center position (x, y) and the angle a of its long axis with respect to the 
laboratory frame. The beating of each flagellum is characterized by a single periodic phase variable, ipz and 
ifR for the left and right flagellum, respectively. The flagellar shapes shown in different colors were tracked 
from high-speed recordings and correspond to a time-difference of 2 ms. This beat pattern was used for all 
computations. B. Binning of tracked flagellar shapes according to shape similarity defines a flagellar phase 
angle as shown on the left. More precisely, we employed a nonlinear dimensionality reduction technique 
as specified in the Supporting Information to represent each tracked planar flagellar shape as a point in an 
abstract shape space. This representation reveals the periodicity of the flagellar beat and supports our de- 
scription of the flagellar beat as a fixed sequence of flagellar shapes parameterized by a single phase variable 
(p. 
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FIG. 2: Hydrodynamic interactions between the two flagella are weak. A. Model Chlamydomonas cell 
represented by an ensemble of 300 spheres used to compute hydrodynamic friction forces at low Reynolds 
numbers. In our calculations, the model cell was assumed to be far from any surfaces. B. Illustration of hy- 
drodynamic interactions between spheres. A single sphere (labeled 1) moving with velocity v\ y > along 
the y-axis will drag fluid alongside and thus exert a total hydrodynamic friction force F\ y = G\i^ y Vi y > 
on the fluid. If a second sphere (labeled 2) is held fixed close to the first one, it will locally slow down this 
fluid flow. The force F2 required to hold the second sphere equals the force exerted by this sphere on the 
fluid; its y-component F2 y = G2\, yy v\ y < defines a friction coefficient G2i, yy that characterizes hydro- 
dynamic interactions between the two spheres. C. Hydrodynamic interactions between different parts of 
the model cell. Analogous to panel B, one defines a matrix Gij :Vy of hydrodynamic friction coefficients for 
the ensemble of 2 • 14 flagellar spheres and the rigid sphere cluster constituting the cell body that together 
represent a Chlamydomonas cell (inset). Each column of the color coded matrix shows the magnitude of 
hydrodynamic friction exerted by a flagellar sphere (or the cell body), if a single sphere or the cell body 
is moved parallel to the long cell body axis. Off-diagonal entries characterize hydrodynamic interactions, 
which are particularly pronounced along a single flagellum (white arrow), or between one flagellum and the 
cell body (central column). Hydrodynamic interactions between the two flagella are very weak and partly 
screened by the cell body. D. Same as in panel C, but for a recovery stroke configuration. There are weak 
hydrodynamic interactions between the proximal segments of the two flagella (white arrow). All friction 
coefficients shown scale with the viscosity of the fluid, which was taken as the viscosity of water at 20° 
Celsius, r\ = 1 pN ms//jm. 
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FIG. 3: A. For synchronized flagellar beating, we compute saltatory forward swimming with positive in- 
stantaneous velocity during effective stroke beating, and a backward motion during the recovery stroke (left 
panel); this behavior is summarized by cartoon drawings to the very right. A typical experimental velocity 
profile of a Chlamydomonas cell in a shallow observation chamber measured during a cycle of synchronized 
beating is shown for comparison in the middle panel. B. Flagellar asynchrony causes cell body yawing, both 
in theory and experiment. Shown is the instantaneous rotation rate a of the cell body in color code as a func- 
tion of the respective phase of the two flagella. For in-phase synchronized flagellar beating (dashed line), 
the cell body does not rotate (green). For out-of-phase flagellar beating, however, we find significant cell 
body rocking (blue: clockwise, red: counter-clockwise). 
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FIG. 4: Flagellar synchronization by cell body yawing. A. For a free swimming cell (top panel), the 
equation of motion 1 predicts a yawing motion of the cell body characterized by a(t) if the two flagella 
are initially out of synchrony (middle panel). The flagellar phase difference 6 is found to decrease with 
time (lower panel, solid line), approximately following an exponential decay ~ exp(— Xt/T) (dotted line), 
where T is the period of the flagellar beat and A defines a dimensionless synchronization strength. Thus, 
in-phase synchronized beating is stable with respect to perturbations. Dots mark the completion of a full 
beat cycle of the left flagellum. B. To mimic experiments where external forces constrain cell motion, we 
simulated the idealized case of a cell that cannot translate, while cell body yawing is constricted by an 
elastic restoring torque Q a = —ka that acts at the cell body center (upper panel). Again, the two flagella 
synchronize (3 rd panel) with a synchronization strength A that can become even larger than in the case of 
a free swimming as shown here for k = 2 • 10 3 pN fim, which is close to the rotational stiffness for which 
the synchronization strength A is maximal (lowest panel). For very large clamping stiffness k, the cell body 
cannot move and the synchronization strength A attenuates to a basal value A « 0.03, which arises solely 
from direct hydrodynamic interactions between the two flagella (arrow). Parameters: 2it/ojq = 30 ms. 
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Theory Experiment 




FIG. 5: Flagellar beating and cell body yawing are coupled in a bidirectional way. Theory (left): In our 
theory, the beat of the left flagellum generates a torque, which, in the absence of the right flagellum, has to 
be counterbalanced by a yawing motion of cell body, see eqn. 4. This effect is quantified by the coupling 
function v(tp) shown in the upper panel, normalized here by po = (p) : the effective stroke (<p/, ~ 0°) of the 
left flagellum causes the cell to yaw clockwise. Conversely, yawing of the cell changes the hydrodynamic 
friction force that opposes the flagellar beat, which, in our theory, speeds up or slows down the beat, see 
eqn. 2. This effect is quantified by the coupling function p(<p) shown in the lower panel: a counter-clockwise 
yawing during the effective stroke of the left flagellum slows down its beat. The coupling of beating and 
yawing allows for flagellar synchronization in a free swimming cell. Experiment (right): By fitting eqns. 2 
and 4 to experimental time series data, we can recover the coupling functions fi(tp) and i/(ip)/po (1 cell, 
n = 5 time series of 0.5 s duration, gray regions denote meanis.e.). 
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VI. SUPPORTING INFORMATION 



A. Image analysis 

High-speed movies were analyzed using custom-made Matlab software (The Math Works Inc., 
Natnick, MA, USA); our image analysis pipeline is illustrated in figure SI: In a first step, estimates 
for position and orientation of the cell body in a movie frame were obtained by a cross-correlation 
analysis using rotated template images. In a second step, these position and orientation estimates 
were refined by tracking the bright phase halo surrounding the cell. The first and second area 
moments of the cell rim provide accurate estimates for the center of the cell body and its long 
orientation axis: While the tracking precision of the first step amounts to < 500nm for the position 
and a few degrees for the orientation, these values are reduced to < 50nm and < 0.5° after the 
second step, respectively. Special care was taken to reduce any potential bias of the flagellar 
phase on the cell body tracking; for example, the cell rim close to the flagellar bases was obtained 
by interpolation instead of direct tracking. The flagellar base is visible as a continuous, parabola- 
shaped curve that connects the proximal ends of the two flagella; tracking of this flagellar base was 
done by a combination of line-scans and local fitting of a Gaussian line model (step 3). Flagella 
were tracked by advancing along their length using exploratory line-scans in a successive manner 
(step 4). Flagellar tracking can be refined by local fitting of a Gaussian line model. A movie 
consisting of thousand frames can be analyzed in an automated manner within 10 hours on a 
standard PC. Movies from red light illumination conditions were of lower quality and required 
manual correction of the automated tracking results for each frame. 

B. Flagellar shape analysis 

We employ a non-linear dimension reduction technique to represent tracked flagellar shapes as 
points in a low-dimensional abstract shape space. In a first step, smoothed tracked flagellar shapes 
corresponding to one cycle of synchronized flagellar beating (shown in figure 1A) were used to 
define the basis of the shape space. Flagellar shapes can be conveniently represented with respect 
to the material frame of the cell using a tangent angle representation [27, 28]. In terms of this 
tangent angle 0(s), the x(s) and y(s) coordinates of the flagellar midline as functions of arclength 
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s along the flagellum can be expressed as 

x(s) = x(0) + / di cos[a + 0(f)] and y(s) = y(0) + / d£ sin[a + 6(£)\. (SI) 
Jo Jo 

Here, a is the orientation angle of the long axis of the cell body (figure 1A), which implies that 
6(£) characterizes flagellar shapes with respect to a material frame of the cell body. By averaging 
the tangent angle profiles 9(s, t) over a full beat cycle, we define a time-averaged flagellar shape 
characterized by a tangent angle 9(s). To characterize variations from this mean flagellar shape, 
we employed a kernel principal component analysis (PCA) [29]. The kernel used to compute the 
Gram matrix D for the kernel PCA must account for the 2n -periodicity of the tangent angle data 
and was taken as Dij = J Q L ds cos[#(s, U) — 9(s, tj)\. The first three shape eigenmodes account for 
97% of the spectrum of D and are shown in figure S2A. The relative contributions to the spectrum 
read 67% (first mode), 18% (second mode), 12% (third mode). While the first mode 6i(s) (blue) 
describes nearly uniform bending of the flagellum, the second mode 9 2 (s) (green) and the third 
mode 9 3 (s) (red) together comprise the components of a traveling bending wave. 

Next, any flagellar shape can be projected onto the shape space spanned by these three shape 
modes: Given a flagellar midline with coordinates x(s) and y(s), we seek the optimal approxi- 
mating shape with coordinates x(s), y(s) whose tangent angle 9(s) is a linear combination of the 
fundamental shape modes 

9(s) = 9(s) + /3Ms) + (3 2 9 2 {s) + p 3 9 3 (s). (S2) 

The coefficients Pi, (3 2 , [3 3 are obtained by a non-linear fit that minimizes the squared Euclidean 
distance J Q L ds\x(s) — x(s) | 2 + \y(s) — y(s) | 2 . This procedure is robust and works even if flagellar 
shapes could only be tracked partially with tracked length V shorter than the total flagellar length 
L. Note that for non-smoothed flagellar shapes, the tangent angle representations can be noisy and 
are thus less suitable for fitting as compared to x, y coordinates. 

A time sequence of tracked flagellar shapes thus results in a point cloud in the shape space 
parametrized by the shape mode coefficients (3\, (3 2 , (3 3 . We fitted a closed curve to the torus- 
like point cloud, see the solid line in figure S2B. This closed curve represents a limit cycle of 
periodic flagellar beating. Each tracked flagellar shape can be assigned the 'closest' point on 
this limit cycle, i.e. the point for which the corresponding flagellar shape has minimal Euclidean 
distance. By choosing a phase angle parametrization for the limit cycle, the phase angle of each 
flagellar shape is determined modulo 2n. A time-series of flagellar shapes thus yields a time-series 
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of the flagellar phase angle (p(t). The phase angle parametrization of the limit cycle had been 
chosen such that the flagellar phase angle tp and its time derivative are not correlated. Finally, the 
zero point y? = was chosen such that the corresponding flagellar shape was nearly straight and 
perpendicular to the long cell axis. 

C. Computation of hydrodynamic friction forces 

For our hydrodynamic computations, we represented a Chlamydomonas cell by an ensemble of 
N = 300 equally sized spheres of radius a = 0.25 pm. The cell body was chosen spheroidal and 
is represented by 272 spheres that are arranged in a symmetric fashion to retain mirror symmetries. 
Each flagellum is represented by a chain of 14 spheres that are aligned along a flagellar midline 
with equidistant spacing. The shapes of the flagellar midlines depend on respective phase angles 
(fL and (fR for the left and right flagellum. These flagellar shapes were taken from experiment 
for one full period of synchronized beating and are shown in figure 1A. We assume that the 272 
spheres constituting the cell body move as a rigid sphere cluster. Each of the flagellar spheres 
represents a cluster with just one sphere, which results in a total of n = 2 • 14 + 1 = 29 sphere 
clusters. We then computed the 6n x 6n grand hydrodynamic friction matrix G for this ensemble 
of n spheres clusters using a freely available hydrodynamic library based on a Cartesian multipole 
expansion technique [23]. Recall that the grand hydrodynamic friction matrix G relates the forces 
and torques exerted by the Qn sphere clusters to their translational and rotational velocities [22] 

Po = G • q . (S3) 

Here, q denotes a 6n-vector that combines the translational and rotational velocity components 
of the n sphere clusters, 

qo = (vix, vi y , v lz , w lx , uj ly , u lz , u nz ), (S4) 

while the 6n-vector P combines the components of the resultant hydrodynamic friction forces 
and torques, 

Po = (F lx ,F ly , Fi z ,T[ x ,T[ y ,T[ z , . . . ,T' nz ). (S5) 

(Primed torques represent torques with respect to the center of the respective sphere cluster.) Fig- 
ure 2C,D in the main text shows a sub-matrix of the grand friction matrix, which was defined as 
Gij :V y = G' 6 i_4 i6 j_4, i, j = 1, . . . , n. In this figure, it was assumed that the long cell body axis 
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is aligned with the y-axis of the laboratory frame, i.e. a — 0, which implies that the sub-matrix 
relates motion in the direction of the long cell axis and the hydrodynamic force components pro- 
jected on this axis. 

For our hydrodynamic computations, the multipole expansion order was chosen as 3. An es- 
timate for the accuracy of our computation could be obtained by increasing the expansion order 
parameter, which changed the computed friction coefficients by less than 1 %. Initial tests con- 
firmed that the friction matrix of only the cell body gave practically the same result as the analytic 
solution for the enveloping spheroid [30]; similarly the computed friction matrix of only a single 
flagellum matched the prediction of resistive force theory [31] assuming a flagellar radius equal to 
the sphere radius. Note that the precise value of the flagellar radius is expected to affect hydrody- 
namic friction coefficients only as a logarithmic correction [5]. 

Below, we consider an extension of the theoretical description given in the main text that addi- 
tionally considers the possibility of an elastically anchored flagellar base, which allows for pivoting 
of the flagellar basal apparatus, see figure S8. In this case, the flagellar midlines were rotated by 
an angle ip. 

A set of 2400 pre-computed configurations was then used to construct a spline-based lookup- 
table of the (reduced) hydrodynamic friction matrix as a function of the degrees of freedom 
and ip. The interpolation error was confirmed to be on the order of 1 % or less. This lookup-table 
was then used for the numerical integration of the (stiff) equations of motion 1 and S17. 

D. Generalized hydrodynamic friction forces 

We employ the framework of Lagrangian mechanics of dissipative systems [32] to define gen- 
eralized hydrodynamic friction forces and derive an equation of motion for the effective degrees 
of freedom in our theoretical description of Chlamydomonas swimming and synchronization. The 
6n degrees of freedom q for the n sphere clusters used in our hydrodynamic computations are 
enslaved by the 5 effective degrees of freedom in our coarse-grained theory, see figure 1. Below, 
one more degree of freedom ip is introduced to characterize pivoting of an elastically anchored 
flagellar basal apparatus. We thus have 

qo = qo(q), (S6) 
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where we introduced the 6-component vector q = (x, y, a,(fL, ¥>r, VO mat comprises the 6 effec- 
tive degrees of freedom. The reduced 6x6 hydrodynamic friction matrix T for these 6 effective 
degrees of freedom can be computed from the grand hydrodynamic friction matrix G as 

r = L T • G • L (S7) 

with a Qn x 6 transformation matrix L given by [12] 

Lij = dqo t i/dqj. (S8) 

The rate of hydrodynamic dissipation can now be equivalently written as a quadratic function of 
either q or q 

K = Qo • G • qo = <f • T • q. (S9) 

The generalized hydrodynamic friction coefficients are depicted in figure S3. In this context, 
generalized hydrodynamic friction forces can be defined as 

Pj = T jx x + T jy ij + T ja a + T jL ip L + r jR <p R + F^ip, j = x, y, a, L, R, ip. (S10) 

Interestingly, the generalized hydrodynamic friction force conjugated to one degree of freedom 
depends also on the rates of the change of the other degrees of freedom, which implies a cou- 
pling between the various degrees of freedom. This fact is illustrated by figure S4. Panel 
A depicts the translational velocities of the flagellar spheres caused by pure yawing of the 
cell body with rate a. This motion is characterized by a 6n-vector of velocity components, 
qo"- 1 = L • (0, 0, a, 0, 0, 0) T . Similarly, the beating of the left flagellum induces hydrodynamic 
friction forces as shown in figure S4B. The resultant force (and torque) components are combined 
in the 6n-vector = G ■ L ■ (0, 0, 0, ip L , 0, 0) T . Figure S4 indicates that the scalar product 
qo"- 1 • Pq L ^ = aT aL ip L does not vanish, which implies a non-zero friction coefficient T aL and thus 
a coupling between cell body yawing and flagellar beating. 

In our theoretical description, the phase dynamics of the left flagellum, say, is governed by a 
balance of the generalized hydrodynamic friction force P L and an active driving force Q L ; sim- 
ilarly, Q R = P L for the right flagellum. In the case of free swimming, force and torque balance 
imply P x = P y = and P a = 0. Together with an equation for P^, these equation allow to self- 
consistently solve for the rate of change q of the 6 degrees of freedom. If one degree of freedom 
were constrained, qj = 0, the corresponding force equation becomes void, since a constraining 
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force Qj equal to Pj then balances the generalized hydrodynamic friction force Pj associated with 
this degree of freedom. 

In general, the active driving forces Q L and Q R will depend on the flagellar phase. This phase- 
dependence is fully determined by the requirement that the flagellar phase speeds should be con- 
stant, (pj = lu , in the case of synchronized flagellar beating with 5 = 0. Here, cu denotes the 
angular frequency of synchronized flagellar beating. Explicitly, we find 

Ql{vl) = u \r LL (^ L ^ L ) + T L R(^ L ^ L )-2Tl y (^ L ^ L )/v yy (^ L ^ L )}. (sii) 

An analogous expression holds for Qr(<Pr). Note that the generalized active driving forces are 
conjugate to an angle, and therefore have the physical unit pN /xm/rad. These phase-dependent 
active driving forces can be written as potential forces Qj = —dU/dtpj, j = L,R, where the 
potential U reads 

/<pl rfn 
dipiQiiVL)- dip R QR(ip R ). (S12) 
•oo J — oo 

The potential U continuously decreases with time, indicating the depletion of an internal energy 
store and the dissipation of energy into the fluid during flagellar swimming. The rate of hydrody- 
namic dissipation equals the rate at which potential energy is dissipated 

K=-U = Q L <p L + Q R <p R . (S13) 



E. Analytic expression for the flagellar synchronization strength 

We present details on the derivation of equations 5 and 6 for the synchronization strength A 
in the case of the reduced equations of motion 2-4. We assume equal intrinsic beat frequencies, 
wl = ujr = cjo and a small initial phase difference, < 5(0) <C 1. To leading order in 5, we find 
relations that link the rotation rate a and the rate 5 at which the phase difference changes, 

ka + p((f, <p)a = -d[v(ip)5]/dt (S14) 
6 = -2fji((p)a. (S15) 

Here (p w cu t denotes the mean flagellar phase. The first equation describes how flagellar asyn- 
chrony causes a yawing motion of the cell body, while the second equation describes how this 
yawing motion then changes the flagellar phase difference. In the absence of any elastic constraint 
for yawing, k = 0, we can solve for 5 

(p - 2fiu)5 = 2fiu'to 5. (S16) 

27 



Now, equation 5 follows from equation S16 using A = — dt5/5 and a variable transformation 

(p(t) =u t + O(5). 

In the case of a very stiff elastic constraint with k 3> pw , we make use of the fact that variations 
of the phase difference 5 during one beat cycle will be small compared to its mean value 5 = (5) . 
As a consequence, equation S 14 can be approximated as ka = —v'uq5q. Using this approximation 
and equation S15, equation 6 follows. 

F. Comparison of experiment and theory 

We can compare instantaneous swimming velocities predicted by our hydrodynamic computa- 
tion with experimental measurements and find favorable agreement, see figures 3 and S5. Note 
that wall effects present in our experiments, but not accounted for by our hydrodynamic compu- 
tations, are expected to reduce translational velocities (but less so rotational velocities) [18]. The 
hydrodynamic computations are based on a fixed flagellar beat pattern parametrized by a flagel- 
lar phase angle, which was obtained experimentally for one beat cycle with synchronized beating 
(shown in figure 1A). The good agreement between theoretical predictions and experimental mea- 
surements for the instantaneous swimming velocities further validate our reductionist description 
of the flagellar shape dynamics by just a single phase variable for each flagellum. Next, we tested 
the applicability of the reduced equations of motion eqs. 2-4 in the experimental situation. For 
this aim, we reconstructed the coupling functions fi((p), v(ip) and p(ip) from experimental time 
series data for a, (p L , and (p R . The coupling functions were represented by truncated Fourier se- 
ries and the unknown Fourier coefficients determined by a linear regression of equation 2, 3, or 
4, respectively, see figure S6. Repeating this fitting procedure for data from 6 different cells gave 
consistent results, see S7. Moreover, the phase-dependence of the fitted coupling functions agrees 
qualitatively with our theoretical predictions. Note that our simple theory does not involve any 
adjustable parameters. 

G. An elastically anchored flagellar basal apparatus 

In the main text, we had assumed for simplicity that the flagellar base is rigidly anchored to 
the cell body. While the proximal segments of the two flagella are tightly mechanical coupled 
with each other by so-called striated fibers to form the flagellar basal apparatus, the flagellar basal 
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apparatus itself is only connected to an array of 16 long microtubules spanning the cell [33]. We 
now consider the possibility that this anchorage allows for some pivoting of the flagellar basal 
apparatus as a whole by an angle ip, see figure S8A. In addition to the five degrees of freedom 
of Chlamydomonas beating and swimming considered in the main text (see figure 1), we now 
include this pivot angle ip as a 6 th degree of freedom. The rate of hydrodynamic dissipation is 
now given by 1Z = xP x + yP y + dP a + 0lPl + <&rPr + i>P^, with P^ being the generalized 
hydrodynamic friction force conjugate to the pivot angle ip. Assuming Hookean behavior for the 
elastic basal anchorage with rotational pivoting stiffness k, we readily arrive at an equation of 
motion that reads in the case of free swimming 

(x, y, a, <p L , tp R , i>) T = T-^O, 0, 0, Q L , Q R , -k^) T . (S17) 

Figure S8B shows flagellar synchronization for a free swimming cell with elastically anchored 
flagellar base: Although some basal pivoting occurs as a result of flagellar asynchrony, the swim- 
ming and synchronization behavior is very similar to the case of a rigidly anchored flagellar base 
as shown in figure 4A. For a cell that can neither translate nor yaw, however, the situation is dif- 
ferent, see figure S8C. We find strong flagellar synchronization provided the elastic stiffness k is 
not too large. Flagellar synchronization by basal pivoting is thus effective also for a fully clamped 
cell. In contrast, for a rigidly anchored flagellar base, the synchronization strength A would be 
relatively weak in this case, being due only to direct hydrodynamic interactions between the two 
flageila. 

Flagellar synchronization by basal pivoting is conceptually very similar to synchronization 
by cell body yawing as discussed in the main text. In the case of a fully clamped cell, we can 
approximate the synchronization dynamics by virtually the same generic equation of motion as 
eqs. 2-4, when we substitute ip for a 

ip L =u -Jl((p L )ip, (S18) 
(pR = uo + ~p(<pR)ij>, (S19) 
kip + p((Pl, ¥r)?P = -H<Pl)<Pl + v(<Pr)<Pr- (S20) 

Here, the coupling functions JL, V and p play a similar role as the previously defined /i, v and p 
for equation 2-4 and show a qualitatively similar dependence on the flagellar phase, see figure 
S9. To derive equations S18-S20, we neglected direct hydrodynamic interactions between the two 
flageila and approximated the active driving forces by Ql{v) = u T LL ((p, (p)\^,=o and Qn{<p) = 
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uoTrr(<P,<P)\1,=o- The coupling functions are defined as pz(ip) = -T^ip, ip)/r LL (ip, v?)^ =0 , 
v(ip) = —T^ L ((p, (p)\^,=o, and p(ip, (p) = T^((p L , (Pr)\^ = q. This choice retains the key nonlinear- 
ities of the full equation of motion, see also figure S3. Equation S18 states that pivoting of the 



speeds up the right flagellum). For synchronized flagellar beating, there will be no pivoting of 
the flagellar base. For asynchronous beating, however, the flagellar base will be rotated out of its 
symmetric rest position by an angle ip if the stiffness k is not too large. Any pivoting motion of 
the flagellar base during the beat cycle changes the hydrodynamic friction forces that oppose the 
flagellar beat, which in turn can either slow down or speed up the respective flagellar beat cycles, 
and thus restore flagellar synchrony. 

To gain further analytical insight, we study the response of the dynamical system SI 8-S20 after 
a small perturbation < 5(0) <C 1. To leading order in 5 = ipi, — <Pr, we find (with ip u t) 



In the biologically relevant case of a relatively stiff basal anchorage of the flagellar basal apparatus 
with k 3> pujQ, we find for the synchronization strength a result analogous to equation 6 



flagellar basal apparatus with ip > slows down the effective stroke of the left flagellum (and 



ktp + p(ip,ip)ip 



d[v((p)S\/dt, 



(S21) 



5 = - 



(S22) 




(S23) 
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VII. SUPPLEMENTARY FIGURES 



0. Raw movie frame 1. Rotated templates 

777777] 




2. Cell rim detection 3. Flagellar base 




4, Flagellar tracking 5. Final result 




FIG. SI: Image analysis pipeline used to automatically track planar cell position and orientation as well 
as flagellar shapes in high-speed movies of swimming Chlamydomonas cells. 0. A typical movie frame 
1. Rotated template images used for a cross-correlation analysis to estimate cell position and orientation 
in a movie frame. 2. The cell body outline was tracked by detecting intensity maxima (green) of line 
scans along rays (shown in blue), which emanate from the putative cell body center. From the cell body 
outline, we obtain refined estimates for cell position and orientation. 3. The position of the flagellar base 
was then determined using a fan of line-scans (along the blue lines), followed by a line-scan (green) in a 
direction perpendicular to the maximal intensity direction (red). 4. Finally, flagellar shapes were tracked in 
a successive manner using similar combinations of line-scans as in step 3. 
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FIG. S2: We represent a single flagellar shape by n = 3 shape coefficients as a point in an abstract shape 
space that is spanned by 3 principal shape modes. A. The principal shape modes were determined by em- 
ploying a kernel principal component analysis (PCA) to the tangent angle representation 9(s) of smoothed 
flagellar shapes that were tracked from the left flagellum of cell no. 2 during one beat cycle of synchronized 
flagellar beating. From the PCA, we obtained three dominant shape modes with respective tangent angle 
representations 9\(s), 9^ (s), 9% (s) as shown. Together, these principal shape modes account for 97% of the 
variance of this tangent angle data set. For sake of illustration, exemplary flagellar shapes corresponding to 
the superposition of the mean flagellar shape and just one shape mode with tangent angle 9{s) + f3i9i(s), 
i = 1, 2, 3 are shown to the right (— 5 < Pi < 5). B. Each tracked flagellar shape from one flagellum can be 
represented by a single point in an abstract shape space that is spanned by the three principal shape modes. 
More specifically, the coordinates (fii, P2, P3) of this point are obtained by approximating the tracked flag- 
ellar shape by a superposition of a previously computed mean flagellar shape and the three principal shape 
modes, see equation S2. The set of flagellar shapes from an entire experimental movie thus corresponds to 
a point cloud. This point cloud scatters around a cl(3&d curve (solid line), which reflects the periodic nature 
of the flagellar beat. This closed curve has been obtained by a simple fit to the point cloud of flagellar 
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FIG. S3: Generalized hydrodynamic friction matrix Fij associated with the effective degrees of freedom 

x, y, a, tfL, ifR, tp. This generalized friction matrix determines the generalized hydrodynamic friction 

forces Pi conjugate to the degrees of freedom q = (x, y, a, tfL, <Pr, VO as Pj = Fijqj, and is computed as 

a projection of the grand hydrodynamic friction matrix, see equation S7. Each friction coefficient Fij is a 

periodic function of the two phase angles (pi and pr, Tij = V^j{ipi,, ipn) and is represented as a colorplot 

with axes as indicated. Here, a is set to zero; different values of a would correspond to a simple rotation of 

the matrix shown. By Onsager symmetry, Tij = Fji [34]. Several features are note- worthy: The coefficient 

T lr characterizes hydrodynamic interactions between the two flagella, and is found to be small compared 

to e.g. Tll. The other coefficients F^j = FjL, which set the friction force Pl conjugate to tpz, depend 

strongly on ipi, but almost not on ipR. This is yet another manifestation of the fact that direct hydrodynamic 

interactions between the two flagella are comparably weak. Analogous statements hold for the coefficients 

Fjij. A counter-clockwise rotation of the cell, a > 0, will increase the friction force Pl during the effective 

stroke of the left flagellum {Fi a > 0), but decrease the corresponding the respective friction force Pr for 

the right flagellum during its effective stroke (FR a < 0). Mirror-symmetry of the swimmer amounts to 

invariance of the friction matrix under the substitution (x, y, a, (pi, (Pr) — > (—x, y, —a, <pr, <pl), which 

implies a number of symmetry relations, e.g. p = F aa must be symmetric in ipL and pr. Finally, this 

rotational friction coefficient p = F aa depends on the flagellar phases in a more pronounced way than 

the translational friction coefficients F xx , F yy . This is inline with the general fact that rotational friction 

coefficients depend stronger (as ~ Z 3 ) on the effective linear dimension / of an object than translational 

friction coefficients (~ I). The coefficients Fj a and Fj^ associated with yawing of the whole cell and 

pivoting of the flagellar apparatus, respectively, show a similar dependence on the flagellar phases. 
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FIG. S4: Coupling of cell body yawing and flagellar beating. A. Translational velocities of the flagellar 
spheres used in our hydrodynamic computation associated with a pure yawing motion of the cell body with 
rate a. B. Hydrodynamic friction forces exerted by the flagellar spheres (as well as by the cell body), if the 
left flagellum advances along its beat cycle with rate tpL- The generalized hydrodynamic friction coefficient 
T a L that couples cell body yawing and beating of the left flagellum can be computed as a scalar product 
between the velocity profile resulting from yawing and the force profile resulting from flagellar beating and 
is found to be non-zero. Parameters: if£ = ojq, a = 0.2uo, 2tt/u)q = 30 ms. 
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FIG. S5: A. Instantaneous swimming velocity in the direction perpendicular to the long cell axis as a 
function of the flagellar phase angles ipL and tpR. For synchronized flagellar beating (dashed line), this 
velocity vanishes in our theory for symmetry reasons (green). If the two flagella are out of synchrony, 
however, significant sidewards motion of the cell is observed, both in theory and experiment. Note that wall 
effects present in the experiments, but not considered in the computations, reduce translational velocities. 
B. Instantaneous swimming velocity in the direction of the long cell axis, again as a function of the flagellar 
phase angles. 
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A Fitting equation 2: <j> L =ra L -u.((|) L )a 




FIG. S6: The reduced equations of motion 2-4 were fitted to experimental time series data for the yawing 
rate a of the cell, as well as the flagellar phase speeds ip^ and ipR. This provided experimental estimates for 
the phase-dependent coupling functions fj,, v, and p. Specifically, we represented each coupling function 
as a truncated Fourier series and determined the unknown Fourier coefficients by a linear regression using 
eqs. 2-4. A. Linear regression of equation 2. Shown to the left in black is the instantaneous flagellar phase 
speed of the left flagellum ipi (smoothed with a span of 15 ms). Shown in red is a reconstructed phase speed 
ujl — pl{^l)^ that depends on the instantaneous cell body yawing rate a, as well as the intrinsic flagellar 
frequency ojl and phase-dependent coupling function \x l for the left flagellum, which were obtained by the 
fit. The coefficient of determination was R 2 = 30%. The estimate for \xl obtained from this fit is shown 
to the right (blue), together with a theoretical prediction (black), see also figure 5A. B. Linear regression 
of equation 4. Shown to the left is the measured instantaneous yawing rate a of the cell body (black), and 
a yawing rate reconstructed from the flagellar phase dynamics, \—v(^p£)^>i, + v{<Pr)<Pt<\/ p(tyLi <-Pr) (red). 
The coefficient of determination was R 2 = 93%. From this fit, we obtain an experimental estimate for the 
coupling functions v(ip) and <pr) (blue); theoretical predictions are shown in black. 
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FIG. S7: Experimental fits for the coupling functions p, v, and p introduced in eqns. 2-4 (blue curves, 
shaded regions indicate meanis.e.) and theoretical predictions (black). The coupling functions p, v, and p 
relate flagellar beating and cell body yawing. Fitting results are shown for 6 different cells illuminated by 
either white or red light as indicated. For each cell, we employed n tits using n non-overlapping time series 
of duration 0.4 — 0.5 s with n as indicated. The blue curves represent the average of the fitted coupling 
functions for the n fits; the averaged coefficient of determination R 2 is stated. 
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A pivoting of flagellar basal apparatus 
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FIG. S8: Theory of flagellar synchronization by basal pivoting. A. We consider the possibility of an elasti- 
cally anchored flagellar basal apparatus (red), which allows for pivoting of the basal apparatus (solid lines) 
by an angle %jj from the its symmetric reference configuration (dashed). B. For a free swimming cell, the 
equation of motion S17 predicts both a yawing motion of the cell characterized by a(t) and a pivoting mo- 
tion of the flagellar base characterized by ip(t), if the two flagella are initially out of synchrony. The flagellar 
phase differences 5 is found to decrease with time (solid line), approximately following an exponential de- 
cay (dotted line). This implies that the in-phase synchronized state is stable with respect to perturbations. 
Each completion of a full beat cycle of the left flagellum is marked by a dot. The synchronization behavior 
in the case of an elastically anchored flagellar basal apparatus is nearly identical to the case of a stiff an- 
chorage as shown in the main text in figure 4. C. For a clamped cell that can neither translate nor rotate, the 
flagellar apparatus can still pivot and will do so if the two flagella are initially out of phase. As in the case of 
a free-swimming cell, the flagellar phase difference 5 decays with time, indicating stable synchronization. 
In the case of a clamped cell, the synchronization strength A strongly depends on the stiffness k of the elastic 
anchorage of the basal flagellar apparatus. Parameters: 2it/oj$ = 30 ms, k = 10 4 pN ^m/rad. 
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FIG. S9: Theoretical coupling functions for the case of a pivoting flagellar base. A. A pivoting motion of the 
flagellar base changes the hydrodynamic friction force associated with flagellar beating and thereby speeds 
up or slows down the flagellar beat cycle in our theory. This effect is quantified by a coupling function p, 
see equation S18. B. Hydrodynamic friction associated with pivoting of the flagellar base (and the attached 
flagella) is characterized by a friction coefficient p, see equation S20. This friction coefficient is maximal 
when the two flagella extend maximally from the cell body during their effective stroke. C. The beat of the 
left flagellum causes pivoting of the flagellar base. This effect is quantified by a coupling function v, see 
eqn. S20. 
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